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The instanton production for out-of-equilibrium scalar fields in a thermal bath of photons is 
computed in the thin-wall limit. Such a process accurately describes the melting and the breaking 
of a pion string via bubble nucleation. In a plasma, the topologically unstable pion string becomes 
metastable so that it can decay through quantum tunneling. In the following, the decay of the false 
vacuum into the true vacuum in the core of the string is used to quantify the quantum decay rate of 
the string into two strings. The cylindrically and spherically symmetric instantons are found. As a 
result of the tunneling, the string core melts, vacuum bubbles propagate along the string at almost 
the speed of light and baryon production could occur via skyrmion production. 



I. INTRODUCTION 



Topological defects are widely observed in condensed 
matter: for instance, string-like defects form when the 
vacuum manifold is such that the first homotopy group 
is non-trivial. In the case of the pion string, the circle of 
minima formed by the charged fields, |7T C | = r\, encloses 
the defect line if finite temperature effects are considered 
[l| 0]. In the same way, non-singular point-like topo- 
logical objects, such as skyrmions, are observed in liquid 
crystals and form when the third homotopy group of the 
vacuum manifold Ai. is non trivial: ir 3 (M) ^ 1@ 0. 
Following these examples, the early Universe should con- 
tain such topological defects due to the rapid phase tran- 
sition it undergoes while cooling down. As a consequence, 
topological defects can play an important role in early 
universe cosmology (see e.g. [H~0] for overviews). On one 
hand, particle physics models which yield defects such as 
domain walls which have problematic and unobserved ef- 
fects can be ruled out. On the other hand, string- like 
topological defects may help explain certain cosmologi- 
cal observations. They could contribute to structure for- 
mation or generate primordial magnetic fields which are 
coherent on cosmological scales @ . 

In the following we focus on a specific period of the 
early Universe, namely the QCD phase transition. A big 
challenge in cosmology is to relate the models we have to 
the physics of the standard model of particles. Here, the 
linear sigma model Lagrangian, which describes pions at 
low energy, provides a starting point for our analysis. In 
this set-up, the pion string is topologically unstable due 
to the fact that the vacuum manifold symmetry group, 
0(4), is too large. In a previous attempt to make the 
pion string stable, we found that the presence of a pho- 
ton bath enhances the stability of the string for a certain 
range of temperatures by effectively reducing the vac- 
uum manifold to a circle [l[ 0. However, PJ indicates 
this string is not stable but rather metastable. As a re- 
sult, the system can move to the lowest energy state, the 



true vacuum, through a process known as bubble nucle- 
ation. The quantum tunneling from the false to the true 
vacuum is studied by finding the instanton solution. For 
the pion string, this mechanism also generates new field 
configurations, skyrmions, via pair production which can 
be interpreted as baryon-anti-baryo n pair p roduction in 
the context of the Skyrme model [9| PJJ ■ 

Instanton production and string breaking have been 
widely studied in literature, mainly from a "world-sheet" 
point of view 

EMI. Here our approach is different and 
more closely resembles more the study of 1st order phase 
transitions in the early universe (l8l-l22j. We are inter- 
ested in the pion string below the chiral phase transition 
temperature of QCD, and thus we only deal with out- 
of-equilibrium scalar fields. Considering a thermal bath 
of photons coupled to the system yields a metastable ef- 
fective potential and allows the presence of embedded 
topological defects (see [22 
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topological defects). Our goal is to compute the tunnel- 
ing probability from the false vacuum to the true vacuum 
of the theory in order to quantify the decay rate of the 
string through bubble nucleations. The string will be 
truly stable and thus have an impact on cosmology, only 
if the decay rate is bigger than the Hubble parameter H: 
H < r where V is the decay rate of the string. In gen- 
eral, different bubble shapes are possible and the total 
decay rate should be the sum of the decay rate into all 
possible configurations. In order to estimate the decay 
rate of the string into two strings, we mainly focus on 
the semi-classical "bounce" computation by Coleman in 
the thin- wall approximation [25| . Making these approxi- 
mations allows us to estimate the euclidean action of the 
instanton for different symmetry configurations. 

In the first subsection, we review the pion string and 
introduce its finite temperature effective potential. Due 
to the shape of the potential, the string is stable above a 
certain temperature, metastable at smaller temperatures 
and classically unstable below To. In the following sub- 
section, we set one ansatz for each pair of fields describing 
the string: (fT()]) and (|17p in order to represent the string 
instanton before (t = — oo) and after tunneling (t = 0). 
We then show the resulting euclidean action. 
Keeping cylindrical symmetry of the equations of mo- 
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tion necessary to get straight strings, implies that we are 
left with three possible configurations: two are cylindri- 
cal and one is spherical. In section IV, we study each of 
them in turn: first, the 0(2) symmetric instanton, then 
the 0(3) one, and finally the spherical one. In general, 
the more symmetry an instanton has, the smaller Eu- 
clidean action it has, so that it contributes more to the 
decay rate. Finally, we compare the value of the differ- 
ent instanton contributions, estimate the decay rate for 
precise values of the linear sigma model parameters and 
conclude on how likely the melting process is. 



II. EFFECTIVE POTENTIAL FOR THE PION 
STRING 



gauge field A^, and four scalar fields out-of-equilibrium: 
two combined into one complex charged field, 7r c and two 
neutral fields forming the pion string combined into the 
complex field, <f>. Above T c h, thermal fluctuations are so 
big that the a field condensate goes to zero and there is 
no string solution. In the case at hand, below the confine- 
ment temperature, non-thermal pions dominate so that 
a zero temperature computation for the effective poten- 
tial (|4]) suffices. The number density of pions being very 
small compared to the one of photons, we can consider 
that photons are indeed in thermal equilibrium while the 
pions are not. In addition, from a cosmological perspec- 
tive, the temperature must also be below the freeze-out 
temperature, Tf, in order for the fields to leave thermal 
equilibrium. 



A. The pion string 

As a toy model for the analytical study of the stabiliza- 
tion of embedded defects by plasma effects we consider 
the chiral limit of the QCD linear sigma model, involving 
the sigma field a and the pion triplet 7? = (it , n , 7r 2 ), 
given by the Lagrangian 



Co = \d^a + l -d^3^ -±(a 2 +T? 2 - V 2 ) 2 , (1) 



Considering temperatures below T c h and Tf give rise 
to a new effective potential at finite temperatures [l|: 
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where rf 1 is the ground state expectation value of a 2 + ir 2 . 
In the following, the potential in ([T]) is denoted by Vq. 

Two of the scalar fields, the a and 7To, are electrically 
neutral, the other two are charged. Introducing the cou- 
pling to electromagnetism, it is convenient to write the 
scalar field sector £ of the resulting Lagrangian in terms 
of the complex scalar fields 



tt + = — -j=(n x + in 2 ) , 
v2 



f 



(2) 



According to the minimal coupling prescription we obtain 
C = ^d^a+^d^d^Tr + D+tt+D^tt- -V , (3) 



where = 



= dp - ieA^ and V 

-2\2 



is a Mexican-hat potential: Vq = j(a + 7? — rfY 
Note that the effective pion-photon interactions appear 
through the covariant derivative in ([3]). 
When one considers a specific configuration where 
a 2 + 7To 2 = rj 2 and 7r + = 7r~ = 0, the vacuum manifold 
S 3 reduces to S 1 and string configurations exist but 
are not topologically stable since iri(S 3 ) = 1. From 
now on we will work with 2 complex scalar fields, one 
combining the charged pions, ir c = ir 1 + iir 2 and the 
other representing the uncharged fields, (f> = a + ittq. 

In order to make the pion string stable let us consider 
a thermal bath of photons in contact with the fields. At 
temperatures below the chiral symmetry breaking tem- 
perature for the linear sigma model, T c % = V2ri, we 
assume that there is only one field in equilibrium, the 



where je is the Euler-Mascheroni constant. For this 
potential, the vacuum manifold forms a circle in phase 
space: \(f>\ 2 = rj 2 and |7r c | = 0. As mentioned earlier, 
such a system admits string solutions: in our case, at 
temperature below the confinement scale, an embedded 
neutral string exists. As shown in Fig[I] above the crit- 
ical temperature, Ti, the potential barrier is infinite in 
the charged pions direction and no tunneling occurs. For 
temperatures between T c and Ti, |7r c | = 0, is the true 
VEV so that the pion string is stable. For both these 
cases, the neutral string is quantum mechanically stabi- 
lized. 

However, for temperatures between Tq and the critical 
temperature T c = ±®M , 1 the system exhibits a 

1st order phase transition due to the negative cubic term 
in the potential. As a result, the zero charged pion con- 
densate becomes metastable: |7r c | is in a false vacuum so 
that we expect quantum tunneling. In other words, this 
quantum tunneling occurs when the charged fields tunnel 
to a non-vanishing expectation value in the true ground 
state, |7r c | = v' where v'is the true VEV (remember that 
there are two charged fields so this is actually a " true vac- 
uum manifold", a circle of radius v'). As a consequence, 
the core of the neutral string melts and the string breaks 
at various locations. Indeed, these quantum tunneling 
effects can happen anywhere on the string, resulting in 
the nucleation of bubbles at the breaking points. Let us 
now take a closer look at the potential giving rise to this 
quantum tunneling. 
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B. Potential in the decoupling limit at 
high-temperatures 



In the high-temperature expansion, 4^ <C 1 so that 
we can neglect higher-order terms of the full effective 
potential (0]) and drop the logarithmic part: 
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Note that for simplicity, we also drop the purely temper- 
ature dependent part of the potential, ■ 

The crucial point is that there are four fields so that 
the tunneling process is, in general, extremely difficult to 
study. Restricting ourselves to the core of the string, we 
can consider = during the tunneling process since 
by continuity, the neutral fields at the center of the string 
vanishes. In its core, the string does not have any poten- 
tial barrier in the <j) direction since <f> vanishes identically. 
It is thus well justified to consider quantum tunneling 
only along the 7r c direction. In other words, it is enough 
to consider the tunneling of only the charged fields, keep- 
ing the neutral ones equal to zero. This greatly simpli- 
fies our task from a computational point of view since, 
as show in FigfTJ the potential becomes a function of ir c 
only. 



where the coefficients are given by 
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and where, for simplicity we set the initial potential en- 
ergy to instead of In order to know when tun- 
neling can occur, we will find at which temperature a 
true vev exists, and above which temperature the barrier 
is infinite. Neglecting Hubble expansion for now, quan- 
tum tunneling occurs at constant energy so that the field 
emerges at an exit point, ir* after which it classically rolls 
down the potential as show in FigJTJ 

In order to find the cxtrcma of the potential we impose 
= and find, in addition to the false vacuum 7T r = 0. 
two other solutions 
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The first point, 7r c +, is a maximum and the sec- 
ond one, 7r c — is a minimum which corresponds to the 
true vacuum. Note that these solutions exist only for 
< 8D (If - T 2 ) , meaning T < T Y with if = 
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In order to get the exit point after tunneling, tt* 
search for the zeros of the potential: 
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Whenever the potential vanishes at 3 points, there will 
be a decay to the lower energy vacuum. This is the case 
whenever 



T > 
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FIG. 1: Finite temperature effective potential in the core of 
the string. The horizontal dashed line indicates the tunneling 
of the string from its initial configuration (where the charged 
fields vanish everywhere) to the exit point after quantum tun- 
neling. The arrow toward the minimum indicates the classical 
decays towards the true vacuum which occurs after quantum 
tunneling. 

In order to study the potential ([5]) it is useful to rewrite 
it in terms of powers of the charged field 
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is the critical temperature. As expected this critical tem- 
perature is smaller than T\. As a consequence, the decay 
into two strings for T > T c won't occur and the initial 
single string configuration remains stable. 
At T = T c , the potential has two degenerate solutions 
7r c = and 7r c = ir® . Moreover, since the two vacua have 



the same energy, V (tt c ) 
the degenerate vacuum: 



0, the exit point, n*, becomes 
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This configuration is the starting point of our instanton 
computation which involves the thin wall approximation. 

Up to this point we have not yet shown which string 
configuration corresponds to the instanton. In attempt 
to do so, we now introduce an ansatz for the pion string 
and show the resulting action. 



III. STRING INSTANTON 




Instantons are classical solutions to the euclidean equa- 
tion of motion with non-zero action. In order to study 
the quantum stability of the string against breaking, one 
can find the instanton configuration which corresponds 
to the tunneling from the false vacuum to the true one. 
In order to do so we first introduce an ansatz for the 
string solution. 



A. Ansatz for the string cylindrical instanton 



FIG. 2: Profile function for the Nielsen-Olesen string for wind- 
ing number n = 1. Contrary to the case for gi(p) and (72 (s), 
the derivative of f(p) does not vanish at the origin p = but 
would do so for n > 1. 



n cb = n for 2 = 0, p < po 

itch = for z ^ and z = 0, p > po 



(16) 



As an example, we present an ansatzt for the cylin- 
drical string instanton which, as indicated by the name, 
has cylindrical symmetry around the string axis. Note 
that the cylindrical symmetry is the minimal symmetry 
required from the string equations of motion. 

Let us consider an initial U(l) neutral string config- 
uration lying on the z-axis which decays to two neutral 
U(l) strings. In our set-tup the initial configuration cor- 
responds to infinite time: t — ±00 whereas the tunneling 
occurs at t = 

t : ±00 -> 

(</>i,7T ci ) y (06,7T c6 ) (13) 

As a result, the solution "bounces" at t = [25| . 

Starting with the Nielsen and Olesen static string [2(| 
at t = ±00 yields 

(<t>i,*ci) = (vf{p)e ine ,0) (14) 

where p is the distance to the core of the string, n the 
winding number and 9 is the angle in the plane orthogo- 
nal to the string. 

After tunneling, the string breaks so that two half- 
strings extend on each side of the z=0-plane. This inter- 
mediate configuration also called, bouncing point, occurs 
at t = 0. As shown in FigJ5J after tunelling, the field </> 
vanishes on the z = 0-plane by continuity in the vicinity 
of the z-axis since it also vanishes in the core of the initial 
and final strings. In order to keep the cylindrical sym- 
metry of the system, we choose the patch where = 
to be a disk of radius po centered on the origin 

fa = Uf(p) einB for z^0andz = 0,p>p Q 
I for z = and p < po 

Thus, for the whole configuration to lie on the vacuum 
manifold around the string, the charged fields have to be 
turned on on the patch defined above 




FIG. 3: Breaking of a string into two strings 

The vertical axis represents the core of the string where the 

field <f> vanishes. In the region of radius po around the core 

of the string of the z = plane, the charged pions turn on 

and take on their v.e.v, (tt c ) — n, while the field <fi is still 

vanishing. 

Let us now introduce a time dependance, assuming an 
extra 0(2) symmetry in the (z,r)-plane where, r = it, 
is the euclidean time. Our ansatz can now be written in 
terms of the new polar coordinate s = \Jt 2 + z' 2 



( t, p, z, 9)= vf(p)e in6 [9i(p) + Vl-9i(p) 2 92(s)\l7) 
tt c (t, p, z, 6)= V y/1 - gi(p) Vl - 92(s) 2 (18) 

Note that two 0(2) symmetries of the bounce solution 
coexist: one in (z, r)-plane and one in the (p, #)-plane so 
that (f)(T, p, z, 9) = (f>rp, s, 6). Rotating back to Minkowski 
spacetime, yields one 0(1,1) symmetry and one 0(2) one. 

The general boundary conditions to get a bounce are: 

9 r (0(r),7r c (r)) T=o = d T {dp b ,Ti cb ) = (0,0) (19) 
lim (&7r c ) = (&,7r cj ) = (0,0) (20) 
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The boundary conditions, (|T? ]) ([T5 ]) ([T9 |) ([20 |) . yield the 
following conditions for gi(p) and g2(s): 

gx (0) = 9i (p < po) = ,gi(p >po) = l ,g[(0) = (21) 
52 (0) = 0, 52 (±oo) = 1,^(0) =0(22) 

where ' denotes the derivative with respect to s (respec- 
tively p) since (|19|) holds at any point of space and we 
can thus substitute r with s (respectively p). 



& (pi g 2 (s) 



7T c = o v ~ c t 
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FIG. 5: Cylindrical shape for the instanton configuration 
When one substracts the bounce action from the original 
string one, one is left with a cylinder of radius po and height 
so that is expanding vertically at the speed of light. 



1 --- 




FIG. 4: Profile function for gi(p) and for g^ (s) 
The profile function for gi(p) and <?2(s) look the same. The 
initially value of gi cancels out and then jumps For the thresh- 
old value p > po, the function goes to I in order to have the 
charged pions turned on only inside the disk of radius po at 
the bouncing point (t=0). The initial value of g2 vanishes 
since (f>(0) = and has a vanishing gradient corresponding to 
the bouncing point as seen in (|22[) . It then jumps to I to come 
back to the single neutral string configuration at t = +oo. 
In a similar way, we asymptotically reach the value I for the 
profile of each of the three functions /, gi and g2 . 



speed of light. As a consequence, the two newly formed 
strings quickly move away from each other along the z- 
axis at the speed of light and the core of the string melts. 

Let us now quantify the probability of the string break- 
ing through quantum tunneling by computing the instan- 
ton euclidean action. 

B. Euclidean action 

Our initial problem includes four fields so that an infi- 
nite number of field space trajectories exist for the quan- 
tum tunneling problem. As a result, we need to choose 
a plausible trajectory. Here we introduce the decou- 
pling limit in more precise terms. Let us define the new 
variable u: u = \f$ 2 ~- j-7rf so that our initial fields are 
4> = usin(6) and tt c — ucos(6). In a schematic way, the 
quantum tunneling description is 



Let us evaluate the typical length scale of the quantum 
tunneling. To go from the initial configuration to the 
bounce one, the fields jump quantum mechanically dur- 
ing a very short transition time that one can estimate 
using Hcisenbcrg uncertainty principle : Et ~ h. The 
typical energy we consider is the energy of the core of 
the string: E ~ jn 4 V where V is the volume. Consider- 
ing the width of the string, w ~ (V^y) -1 , the volume of 
the initial cylinder inside the string is, V ~ irw 2 so where 
so is the height of the cylinder. We therefore obtain the 
typical tunneling length, I ~ t ~ £rf^ a m units where 
H = c = 1. For so big enough, I < w and the assump- 
tion that tunneling occurs within the core of the string 
is consistent. Here so corresponds to the point where the 
quantum jump occurs, meaning that we can infer that 
the profile of 172 (s) has a sharp change at s = so. In a 
similar manner, we suppose that the profile of gi(p) has 
a sharp change at p = po as show in Fig 2] Later, in 
section IV. A, we find the exact relationship between p 
and so in the thin-wall approximation case. 

The effects of this quantum tunneling is to split the 
string. Indeed, the boundary of the cylindrical region 
separating the two strings expand vertically in a way that 
z 2 — t 2 = Sq. As a result the velocity of the wall is v = 

/ z 2_ s 2 

^| = — — - which, after a very short time, becomes the 



t : -> 00 

configuration : bounce — > stable string 

(f) = usin{6) : -> nf{p)e inu > 

tt c = ucos{6) : 7T* -> (23) 

so that tanO = — : — > 00. 

The euclidean action for the string is now 

S E = 47 r 2 I Pdpsds\[{ff + (ff 

W[(|) 2 + (|) 2 ] + 2V>,0,T)]. (24) 

where we have neglected the winding gradient of the 
string. We can choose a specific trajectory for 8(p, s) 
that satisfies the boundary conditions ([23]) . Restricting 
ourselves to the core of the string, (f> = at any time. As 
a result, tanO always vanishes in the decoupling limit. 

We can now write the Euclidean action for the cylin- 
drical string instanton described above: 

S E = 4tt 2 / pdpsdsh&f + (^) 2 + - 2 n 2 <?\ + 
J 2 d s d p p 2 

4n 2 J p d psds[\[{^) 2 + (^) 2 ] + V(<f>,« c ,T)] (25) 
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where the integral runs from to oo. According to Der- 
rick's theorem, the energy configuration resulting from 
the gradient term due to the winding around the string, 
■i?n 2 (j) 2 is infinite in more than one dimension. Fortu- 
nately, the introduction of gauge fields cure this problem 
and from this point of view, the action should be free of 
divergence. 

Let us now estimate the analytical value of the instan- 
ton action and the corresponding string decay rate. In 
order to do so, we can use the thin-wall approximation 
close to the critical temperature, T c . 



IV. THIN- WALL APPROXIMATION 

The thin- wall approximation allows to use a semi clas- 
sical computation to estimate the production rate of thin- 
wall bubble and is valid when the potential is almost 
degenerate. For our potential ((6]), the approximation is 
valid for temperatures just below the critical temperature 
T c . Therefore the first instantons formed correspond to 
thin-wall bubbles. At lower temperature, the thin wall 
approximation breaks and the wall becomes thick since 
the height of the barrier is small compared to the value 
of the potential at its minimum. Whether bubble nu- 
cleation occurs via thick or thin wall tunneling depends 
on how big the decay rate is in the thin wall case. If 
the decay rate is big enough, nucleation will proceed via 
thin- wall tunneling. If not, one can solve the thick- wall 
case which most likely requires numerics. 

Let us consider our potential at temperatures where 
the difference in energy between the two minima is neg- 
ligible. Since the 2 minima are almost degenerate, one 
can neglect the damping term in the equation of motion 
for the charged fields, ir c , so that 



d n r 1 dn 



+ irp = V '(tt c ) 



dp 2 r dp 
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V'(tt c ) ~ V£(7T C ) 



(26) 



(27) 



where Vd(ttc) is the potential in the limit where the po- 
tential has an exact degeneracy (e ~ or equivalently, 
T ~ T c ). 

Let us now rewrite the potential (|6]) to make appear 
the thin- wall parameter e, as in [27| 



-7T C (7T C - 7T C ) - -€TT C TT C (28) 



where e 
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Note as a consistency check that this potential vanishes 
at tt®(T c ) since e vanishes at the critical temperature and 



that the critical value tt®(T c ) matches the one found in 

mi. 

For this potential (f28|) , the one-dimensional Euclidean 
action is 



*t[i(^) 2 +n«] 



Making the temperature dependence explicit yields 



6V2 
(31) 



Si CO = 
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[D(T 2 - T 2 )] 3/2 = 



18V6A 



[(T 2 - T 2 )f' 2 
(32) 

Note that this action has mass dimension 3 since it is the 
one-dimensional action. 

Another quantity, the potential energy density differ- 
ence between the two minima, plays an important role in 
the computation of the action, 
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D\A 



AV = V(n cmm ) = -e(TT^) 
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AV = 8(ET - VXDJT 2 - T 2 )[—(T - Tq )] 3 
* A 

Here, contrary to the Mexican hat potential case (with 
no cubic term ) (2f|, AV, is not equal to the thin- wall 
parameter e which is a dimcnsionless quantity that de- 
creases with temperature (see FigJBJ) 



T= 



and e(T> 
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3tt 2 A 
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FIG. 6: Plot of the thin parameter, e, as a function of tem- 
perature for A = 10e 4 , e 2 = 77 = 66A7MeV. These 
parameters correspond to To = 155.936 MeV, and the criti- 
cal temperature equals to T c = 156.2MeV. The thin-wall ap- 
proximation is valid only in a narrow energy range (less than 
lMeV) very close to T c . The horizontal line represents a 25% 
threshold value (found in (|48[) ) above which the thin wall ap- 
proximation is not valid anymore. For the approximation, to 
be valid, e must be much smaller than this threshold value. 
We see that the thin-wall parameter vanishes at the critical 
temperature. 
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From (|29|) and (|34)) . we can see that the thin- wall 
parameter exactly vanishes at the critical temperature. 
Above, T c , the second minimum of the potential is higher 
than the initial one so that the field is already in the 
true vacuum and no tunneling occurs: The thin-wall ap- 
proximation does not make sense anymore and e becomes 
negative. Far below T c , the value of the thin wall param- 
eter explodes and the approximation breaks down. As 
a consequence, the thin- wall approximation in only valid 
in the vicinity of T c . 

In order to evaluate the thickness of the bubble wall, 
one can solve the differential equation (|27|) where the 
damping term has been neglected. The solution close to 
r ~ R (the radius of the bubble wall which we compute 
in the following) is |27| : 



7r c (r) = ^[l-tanh[ 7 ^} 
where <5, the wall thickness, is equal to 




2V6 



T 2 



(35) 



(36) 



and, as expected, the higher the temperature, the thincr 
the wall is. 

Let us now use the thin-wall approximation and the 
physical quantities mentioned above to estimate the ac- 
tion for different instanton configurations. 



The set of (po,so) that fulfills the above extrcmization 
constraint of the euclidean action form a curve in the 
(p, s) - plane: p s = jgrp- 
Plugging it back into ([55| yields 



S, 



cyl 



64tt 



6 2 S 2 _ 64tt 2 e 4 (T 



T Q 2 ) 2 



AV 2781 A 2 



AV 



128tt 2 
9eA 



(39) 



The string will break at some point on the z=0 plane 
with a probability per unit volume 



r 



sphere 

V 



p 



26 



(40) 



where P2 is a prefactor that wc will estimate later. It is 
worth noting that in the thin- wall approximation e«l, 
so that most of the contribution to the decay rate comes 
from the exponential part. 

After, tunneling, the two strings immediately move away 
from each other along the z-direction at almost the speed 
of light. Moreover, the damping effect coming from the 
plasma pressure should not have a great effect since the 
string expand vertically and not radially. 

There is another possibility for a cylindrical instanton: 
one that would expand radially rather than vertically (see 
FigEJ) , the 0(3) symmetric instanton. 



B. Cylindrical Bounce: 0(3) symmetry 



A. Cylindrical Bounce: 0(2)x 0(2) symmetry 

Let us come back to the string cylindrical instanton 
and compute the most probable configuration by using 
the thin-wall approximation described above and by ex- 
tremizing the action. For this cylindrical case illustrated 
in Fig. [SJ the action ([23]) simplifies to 



S, 



cyl 



= 4tt" 



pdps ds[ \[(^r + ( di 

n 2 s 2 
Po s 



*° ) 2 ] + V(n c )} 



^[p s S 2 -(f ■£ AV] 



(37) 



where S2 is the 2-dimensional action that we approxi- 
mated as: S2 ~ 4(551, so that the action becomes 



S C yl = 47r 2 p so[4^S'i 



Po-sq 



AV] 



(38) 



Let us now introduce the 0(3) cylindrical instanton 
which, when formed, would expanding radially in the 
plane orthogonal to the string. In this case, the gradient 
energy from the winding around the string and the pres- 
sure of the plasma should prevent this radial expansion to 
occur as shown in FigJTJ If no such physical effects would 
prevent radial expansion, the 0(3) instanton should give 
a larger decay rate than the 0(2) one because of its larger 
symmetry group. 



\ 

\ 

N 



7Tc = 



v~c 



3 



The first term is the gradient energy associated with the 
cylindrical wall and the second one corresponds to the 
vacuum energy in the interior of the cylinder. 

In order to find the optimal dimension of the most 
likely nucleated cylinders, we look for the saddle points 
of the euclidean action 



FIG. 7: Cylindrical shape of the 0(3) instanton configuration. 
When one removes the bounce action to the original string 
one, one is left with a cylinder of radius po and height so that 
should be expanding radially at the speed of light. The dashed 
walls on each side of the string represents the gradient energy 
from the string winding that prevents the radial expansion. 



Up 



and 



d s 







For the sake of generality, we now present some results 
that go beyond the thin- wall approximation [22|. 
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The cuclidcan action for an 0(3) bounce for purely ther- 
mal fields reads 



Y = 4 - 85 -^r/(«) 



(41) 



1 + 



2.4 



0.26 



where f(a) = 1 + § 

M 2 (T) = 2D(T 2 -T 2 ) is the mass term and a = ^jffffl ■ 
The approximation above is valid for any temperature be- 
tween To and T c so that f(a) is correct with an accuracy 
about 2% when a goes from to 1 . Applying (|4*T|) to our 
case yields the dimensionless euclidean action 



5, 



cylrad ' 



4.85 



z M 3 (T) 



/(<*) 



(42) 



E 2 T 2 

where zq is the invariant height of the cylinder nucleated. 

The physics in [22| differs from our pion string case 
where only the photons are in thermal equilibrium. How- 
ever the result for the three dimensional Euclidean action 
still holds since the Lagrangians of both systems have the 
same form. The only difference is that, the 0(3) sym- 
metry applies to the three spatial directions (rotational 
invariance of the bubble nucleated), instead of two spa- 
tial plus one time component for the radially expanding 
cylinder (cylindrical configuration nucleated) . 

In general, the most probable configurations nucleated 
should be bubbles with an 0(4) symmetry in the Eu- 
clidean plane. So, let us now have a closer look at the 
0(4) instanton. 



C. Spherical Bounce: 0(4) symmetry 

In general, the most probable way for the charged fields 
to tunnel through the potential barrier is through the 
nucleation of bubbles (which have an 0(4) symmetry in 
the Euclidean plane) due to the fact that 0(4) is the 
largest possible symmetry group [25[. Here, the spheri- 
cal bubble, which incorporates the cylindrical symmetry 
needed to get strings, should give the biggest decay rate. 
However, the two antagonist effects to bubble nucleation 
mentioned above, plasma pressure and string winding 
gradient, should prevent or restrict any radial expansion. 
Thus, we expect that once the bubbles form, they get 
deformed and propagate along the string. 

Let us now take a closer look at the dynamics of the 
problem. The equation of motion for a spherical bounce 
in 4-dimcnsion 



d 2 ir c 
dr 2 



r dr 



(43) 



becomes the same as the cylindrical case since we do not 
consider the damping term in the thin wall approxima- 
tion: = Vd(tt c ) where r and p are the radial dis- 
tances in spherical and polar coordinates, respectively. 
Integrating over a 3-sphere of radius R yields 



S sph ere = ^ J rHr[\{^) 2 + V(w c )} 

= -4R A AV + 2tt 2 R 3 S 1 



(44) 




FIG. 8: Spherical shape for the instanton configuration 
In the ideal sphere of radius R is generated and ex- 

pands radially at almost the speed of light 



where R is the radius of the bubble in Euclidean space. 
As usual, we should extremize the action in order to de- 
termine the radius of the bubble nucleated: = and 

d R 




Plugging back in the action yields 

St 



"^sphere 



21tx 1 



2AV J 



48Ae 3 



(45) 



(46) 



Replacing e by a function of temperature 
see that ^ is positive for any temperature below T c so 
that the radius of the bubble increases with temperature. 
As a result, when the Universe expands and cools down, 
the pion string would generate bubbles of smaller size. 
Using the action we can finally write the decay rate per 
unit volume: 



sphere 

V 



P4 exp[— 



1 



48 A(^ 



e 4 T 2 
3ir 2 A T 2 -T 2 



(47) 



where P4, the prefactor can be estimated on dimensional 
grounds. Therefore the decay rate decreases with tem- 
perature. 

After tunneling in vacuum, the bubble radially ex- 
pands at the velocity v = = ^^gi which rapidly 
becomes the speed of light. However, in our case, the 
pressure the plasma exerts on the wall slows down this 
expansion. 

Let us write now use the velocity we just found and 
write the energy of the bubble wall in the core of the 
string (otherwise we have to consider the neutral pion as 
well) 

E waU =ATr\x\ 2 (S^)(l-v 2 )-i 



which finally reduces to 

, 1 



En 



ill 



4ir\x[ 



27re, 



As expected, when the temperature decreases, the energy 
of the wall decreases while the initial velocity of the bub- 
ble wall increases due to a smaller initial radius. In ad- 
dition, the pressure exerted by the thermal bath is lower 
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at smaller temperature and thus has a milder effects on 
slowing down radial expansion. 

Let us now discuss our results and compare the 0(4) 
instanton we just described with the 0(2) instanton. 



V. DISCUSSION AND NUMERICAL RESULTS 

A. Prefactor estimate 

In order to estimate the decay rate of the string, we 
made use of the semiclassical approximation which only 
applies when the instanton action is very small, Se 3> 1. 
In addition, this approximation holds only if the radius 
of the bubble wall is much bigger than its thickness 



3Si /8 1 



(48) 



A simple rescaling of variables allows our effective La- 
grangian to resemble the one of the linear sigma model by 
absorbing the electromagnetic coupling constant e 2 into 
a new charged field, tt c , and a new coupling constant A. 
Our effective potential (J6j) then becomes 



Veff(ir c ,T) 



12 ^ 

= (T 2 



6\tj 



6A?7 2 



Ttt 3 



6tt 
6tt 



4 "c 



4 "c 



where n c = en c and A = -4 are the new rescaled vari- 
ables. Note that at zero temperature, this potential is 
exactly like the linear sigma model potential. As a re- 
sult, we can use some typical values for the parameters 
in the linear sigma model: A = 10e 4 , e 2 = j^, ?] = 66.47 
McV. This choice of parameters corresponds to a critical 
temperature of T c = 156.2 MeV and a minimum temper- 
ature of T = 155.936 McV. 



or equivalcntly, e« j or T ^> 



In this regime, the tunneling rate is determined mainly by 
the action, and the value of the prefactor, P, can simply 
be estimated using dimensional analysis. According to 
[IH, the prefactor, P in front of the exponential factor 
for the decay rate is of order 



P„ 




on 



where det' denote the determinant computed with the 
zero eigenvalue omitted and n corresponds to the num- 
ber of zero modes associated with translational invari- 
ance in euclidean space [H)]. For the spherical case, 
n = 4, the cylindrical 0(3) case, n = 3 and the cylin- 
drical 0(2) x 0{2) case, n = 2. S E is the corresponding 
euclidean action for the 0(n) symmetric configuration. 
The prefactor computation for different symmetries of 
the instanton is more precisely discussed in [29[ [3(| . 

In order to evaluate the potential, we compute the 
value of the potential close to the critical temperature 

which, is of order T so 



l—D\ ?VT 2 -Tg 

that the total decay rate per unit volume has unit of T . 
Then, a quick estimate for the cylindrical and spherical 
instanton prefactor (|4l?|) yields 



3.7e' 



r 2 -T 2 anrfP 4 ~ L8 A y e \ / 



Let us now turn to the numerical evaluation of the 
exponential part of the decay rate. 

B. Numerical results 

In order to evaluate the decay rate of the string, let 
us plug in the values of the linear sigma model that give 
rise to the observed pion mass and possible sigma mass. 



TABLE I: Decay rate for A ~ 10 for the spherical (0(4)) and 
cylindrical (0(2)) instantons for two different values of the 
thin-wall parameter, e. 



Symmetry 


Action Se 


r 
V 


for e 


~ 0.1% 


r 
V 


for 


e~3% 


0(2) x 0(2) 


12»tv' a 
9eA 




- P 2 e" 


1.7.10° 






-5.5.10* 


0(4) 


7V 2 

48AE 3 




- P 4 e" 


2.4. 10 9 




Pie 


-9.0.10 4 



As show in Tabic HI the exponential part of the de- 
cay rate has almost the same order of magnitude for the 
spherical and cylindrical case when e is big enough (or 
equivalently, the temperature of the thermal bath of pho- 
tons is closer to To). However, we should not forget that e 
must be small for the thin-wall approximation to remain 
valid. 



C. Discussion 

Let us now compare the decay rates for the spherical 
and the cylindrical case. The ratio of the decay rate is 



sphere 



e Ac L 48Ae 



(51) 



Recalling that the thin wall parameter e is very small, the 
decay rate for both configurations is very small due to the 
strong exponential suppression coming from the - factor. 
Moreover, we see that the cylindrical instanton is more 

likely than the spherical one whenever e < ^\/f ~ 4%. 
This is in perfect agreement with the numerical results 
of Table |IJ However, this is not what is usually expected 
from the symmetry argument and the spherical action, 
which has a bigger action, yields a smaller decay rate 
than the less symmetric 0(2) instanton. 

Let us now consider the cosmological perspective and 
compare the decay rate of the string per Hubble vol- 
ume to the Hubble parameter. In general, if T ~ H 
then a computation at finite temperature is necessary and 
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for r <C H a zero temperature computation is enough. 
Here, the requirement for our string to be stable un- 
der breaking is that y <^ H 4 . Roughly speaking we 
know that H ~ ^— i— — where tQCD ~ 10 _5 s is the typical 
time at which the QCD phase transition occurs so that 
H 4 ~ 10" 20 s- 4 - e' 46 s- 4 . As obvious from TablcH the 
decay rate of the pion strings per unit volume is much 
smaller than H 4 and therefore the string is stable close 
to the critical temperature. 

However, the thin- wall approximation is valid in a very 
narrow range close to the critical temperature, T c . As one 
can see in Fig. [TJ close to T c , the height of the barrier, 
V(tt c +), is large compared to its depth, V(ir c — ) which 
actually vanishes at T = T c . In that case, the field value 
changes very quickly from inside the bubble to outside 
and the wall is thin: -c 1. 

But at lower temperatures, closer to T , the height 
of the barrier of the potential becomes smaller than its 
depth, so that the wall becomes thicker. In that case, the 
height of the barrier is of order the depth of the potential 
and one of the minima is much lower than the other. As 
a result, the total potential energy contained inside the 
bubble is smaller and the surface term dominates since 
the wall is thick. We also expect a smaller action coming 
from the thick wall bubble yielding a smaller argument 
in the exponential part of the decay rate. As a conse- 
quence, for the thick wall, the prcfactor includes most of 
the relevant information for the decay rate and we still 
expect a bigger decay rate than for thin- wall bubbles. So, 
at small enough temperature, the larger rate of bubble 
nucleation can make the string unstable. If it is the case, 
skyrmion configurations can form and expand along the 
string at the speed of light. This mechanism will be fur- 
ther developed in a follow-up paper. In this scenario, our 
string, which is stable at high temperature close to T c , 
becomes less stable at lower temperature to end up with 
a classical decay below Jo- 
in contrast, it is worth noting that in the sphaleron 
case, where only thermal fluctuation effects are consid- 
ered, the colder the system is, the less likely it is to hop 
over the potential barrier, thus making the system more 

and more stable. Indeed, when all the fields are in ther- 

s 3 

mal equilibrium, the decay rate goes like r ~ e t where 
53 is the 0(3) symmetric action so that a decreasing tem- 
perature exponentially suppresses bubble nucleation. 

Finally, in the limit where the barrier height ap- 
proaches zero, close to To, the bubble size approaches 
zero and the field classically rolls down the potential 
without any tunneling required so that our string van- 
ishes. 



VI. CONCLUSION 

We proved that considering plasma effects on the 
charged fields helps stabilize the single string configu- 



ration against its decay into two strings for temperatures 
close to the critical temperature, T c , and so contributes 
to cosmology during a certain time. 

However, when the wall becomes thicker at lower 
temperatures, the resulting instanton action should be 
smaller, making the system less stable. In that case, nu- 
merical computations can be carried out in order to pre- 
cisely compute the pion string decay rate without any 
thin-wall approximation. 

Two main effects prevent the string from decaying into 
radially expanding configurations: Plasma pressure, and, 
more importantly, the string winding gradient which cre- 
ates an energy barrier that encloses the core of the string. 
Neglecting those effects we find that the radius of the 
bubbles nucleated increases with temperature. However, 
after nucleation, we expect the bubble to be deformed 
and to propagate along the string at almost the speed 
of light. As a result, the string core melts and the two 
newly formed string segments move away from each other 
at the same speed. This melting is more likely to occur 
far below the critical temperature, when the string decay 
rate is large enough compared to the Hubble parameter. 

In addition, the decay of the string has some interest- 
ing consequences: Opposite winding direction configura- 
tions for the charged fields components form skyrmion 
anti-skryrmion pairs when they meet. If the two configu- 
rations have the same alignment, they just pass through 
each other and no skyrmions form. When skyrmions 
form, they can spread outside the strings since cosmic 
strings are moving. Interpreting those skyrmions as 
baryons § [l(| [HI > the decay of the pion string can con- 
tribute to the production of baryon-anti baryons pairs. 

Finally, it is a great challenge to compute the full tun- 
neling probability in the case where we consider the four 
scalar fields involved in the linear sigma model. Here we 
restricted ourselves to the core of the string where the 
neutral complex field vanishes but there are other trajec- 
tories in field space to explore in order to quantify the 
stability of the pion string. 
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